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ABSTRACT 

Recent X-ray observations have shown evidence for exceptionally broad and skewed iron Ka 
emission lines from several accreting black hole systems. The lines are assumed to be due 
to fluorescence of the accretion disk illuminated by a surrounding corona and require a steep 
J> ' emissivity profile increasing in to the innermost radius. This appears to question both standard 

CO . accretion disc theory and the zero torque assumption for the inner boundary condition, both of 

which predict a much less extreme profile. Instead it argues that a torque may be present due 
. to magnetic coupling with matter in the plunging region or even to the spinning black hole 

f^i ' itself. Discussion so far has centered on the torque acting on the disc. However the crucial 

CO • determinant of the iron line profile is the radial variation of the power radiated in the corona. 

Here we study the effects of different inner boundary conditions on the coronal emissivity 
. and on the profiles of the observable Fe Ka lines. We argue that in the extreme case where 

Oh a prominent high redshift component of the iron line is detected, requiring a steep emissivity 

q . profile in the innermost part and a flatter one outside, energy from the gas plunging into the 

5— i ' black hole is being fed directly to the corona. 
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U 1 INTRODUCTION 

The much improved sensitivity of the X-ray satellite XMM Newton is quickly expanding our capability of performing accurate X-ray spec- 
troscopy of accretion powered systems. In several recent observations of accreting black holes, extremely broadened and redshifted iron Ka 
lines have been detected. Wilms et al. (2001) presented the observations of the nearby Seyfert 1 galaxy MCG-6-30-15 in a low-luminosity 
state and deduced from the iron Ka line profile that the reprocessing material must be located in the very innermost part of the accretion flow. 
Similarly, Miller et al. (2002), found a very broad, skewed Fe Ka emission line in the galactic black hole XTE 11650-500 in its very high 
state. In both cases, not only is the large observed redshift considered as a signature of a rapidly rotating black hole (so that the innermost 
stable orbit is close to the event horizon and the accretion disc can penetrate much deeper in the potential well), but also the disc emissivity 
needs to be strongly peaked towards the inner boundary in order to fit the observed line profile. For simplicity, the disc emissivity is usually 
assumed to have a power-law profile F(r) oc r- _/3 (where r = R/Rs = Rc 2 /2GM is the distance from the center in units of Schwarzschild 
radii), and the index of such power-law is generally assumed to be given by the accretion disc emissivity index in a standard geometrically 
thin and optically thick disc l Sh akura & Sunyaev 1973) . In this case we expect /3 < 3. The fact that both Wilms et al. (2001) and Miller et al. 
(2002) require a larger value of f3 to fit the line profile has prompted speculation that, in the inner disc region, additional energy dissipation 
must be taking place l Agol & Rrolik 2000 Li 2002 1. This might be due to magnetic torque exerted by the accreting material in the plunging 
region at the inner disc boundary iGammie 1999 Rrolik 1999 1 or to magnetic extraction of the spin energy of the black hole tLi 2002> 
through magnetic field lines threading the ergosphere and the accretion disc. Such findings have been confirmed by the analysis of a long 
observation of MCG-6-30-15 with XMM-Newton reported by Fabian et al. (2002). These authors found that "the strong, skewed iron line is 
clearly detected and is well characterized by a steep emissivity profile within 6GM/c 2 and a flatter profile beyond." 

Here we first make the point that, in the framework of standard, relativistic accretion disc models with zero torque at the inner boundary, 
the emissivity profile of the hard illuminating X-rays is more centrally peaked than that of the usually assumed standard accretion disc. This 
is a general property of the energy deposition rate in a magnetic corona above standard accretion discs. Furthermore, the quality of the X-ray 
spectroscopic data is such that, using our coronal models, we can start considering the nature of the inner boundary condition for relativistic 
accretion flows. Using a simple analytic formulation, first introduced by Agol & Krolik (2000) to include magnetic stresses at the marginally 
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Figure 1. On the left (panel a), coronal emissivity profiles F c (solid lines) and corresponding standard disc emissivities Q (dashed lines), both normalized to 
their maximum level, as a function of the disc radius, in units of the inner accretion disc radius (fixed at the marginally stable orbit) for three different values 
of the dimensionless black hole angular momentum: a = 0; a = 0.9 and a = 0.99, from bottom to top. On the right (panel b), the upper panel shows the 
location of the innermost stable orbit (dotted line), the radius of maximal coronal emissivity (solid line) and that of maximal disc emissivity (dashed line), all 
in units of Schwarzschild radii, as functions of the dimensionless black hole spin parameter a, for the case of no torque at the inner boundary. In the lower 
panel the corresponding radial velocity of the gas (in units of the local sound speed) and the relative departure of the gas angular velocity from the Keplerian 
value (assumed in the computations), at the radius of maximal coronal emissivity are shown. 



stable orbit in standard disc equations, we can estimate the effect of a modified inner boundary condition on the observed iron line profiles. 
In particular, we consider the three following possibilities: 

1) There is no net torque at the inner boundary (no-torque case; NT); 

2) Magnetic field lines couple the plunging region to, and exert stress on, the cold, geometrically thin disc only (disc-torque; DT); 

3) Magnetic field lines couple the plunging region to the vertically extended, magnetically dominated corona only (coronal-thread; CT). 

We show that, while in the NT case the coronal emission is more centrally peaked than the disc one (due to the radial dependence of the 
fraction of energy dissipated there), the effect is greatly enhanced in the CT case, with the strongest deviation from the standard emissivity 
profile appearing in the innermost coronal region. Although the emissivity profiles again become more centrally peaked when the torque is 
applied directly to the thin disc (DT), the coronal power is much reduced, and we expect the emission to be completely dominated by the 
quasi-thermal radiation coming from the optically thick disc. 

The structure of the paper is the following: in section 2 we solve the fully general relativistic structure equations for standard geomet- 
rically thin and optically thick accretion discs, coupled to a magnetically dominated corona (with the coupling discussed in Merloni, 2003, 
hereafter M03) and opportunely modified by the inclusion of a term describing the additional dissipation from magnetic stresses at the inner 
boundary either in the disc or in the corona. In such a way we are able to compute the radial profiles of the hard coronal radiation that illumi- 
nates the cold disc and is responsible for the fluorescent iron Ka line production (section|3J. Then, by means of a simple broken power-law 
approximation for the coronal emissivity, in section|4] we simulate and compare the expected line profiles for the different cases. In section|5| 
we discuss a number of possible shortcomings of the model we have adopted, and in section|5|we draw our conclusions. 



2 CENTRALLY PEAKED CORONAL EMISSIVITY: THE NO-TORQUE CASE 

In Merloni (2003; see also Merloni & Fabian 2002), we have shown how it is possible to build a physically self-consistent analytical model 
for a thin accretion disc-corona system by solving algebraic equations for the fraction of power released into the corona as a function of the 
distance from the central (non-rotating) black hole. Here we extend such a treatment to the general relativistic case of accretion discs around 
rotating black holes by calculating the disc-corona structure equations (presented in the Appendix) which are a generalization of the classical 
ones first derived by Novikov & Thorne (1973), and take into account both the effects of the presence of a corona on the disc structure and 
those of the different viscosity prescription on the coronal emissivity profiles. 

The nature of viscosity in accretion flows is still one of the most fundamental open questions of black hole astrophysics. At present, 
magneto-rotational instability (Balbus & Hawley 1991 Ba lbus & Hawley 19 98 ) is favoured as the primary source of the turbulent viscosity 
needed to explain the luminosities of accreting black holes. For this reason we expect that the turbulent magnetic stresses responsible for the 
angular momentum transport in the disc are given by 



Iron line profiles and inner boundary condition in accretion flows 3 



P W R , ^ < B «* B *> - P mag . (1) 

Once the relationship between magnetic and disc pressure (given either by the gas or by the radiation) is established, the accretion disc 
structure can be fully described. This is a long standing issue in accretion disc theory and, in the framework of MRI-driven turbulent 
viscosity models, can be translated into the uncertainty about the mechanisms by, and the level at, which the disc magnetic field saturates. 
Numerical simulations are crucial to test and verify different hypothesis on the outcome of the instability and the final scaling of the magnetic 
pressure. However, full three dimensional MHD simulations of accretion flows that include radiative transfer and is capable to treat radiation- 
pressure-dominated solutions (on the line of those presented in Turner, Stone & Sano, 2002) are needed, in order to answer the question 
of the viscosity scaling for luminous thin discs. On the other hand, by means of a more treatable analytic approach, Merloni (2003) has 
shown that, as the growth rate of the MRI in radiation-pressure-dominated discs is reduced by a factor c g /vA (Blaes & Socrates 2001 
ITurner, Stone & Sano 2002 1, where c g is the gas sound speed and va is the Alfven speed, and if the main saturation mechanism for the 
magnetic field is vertical buoyancy (as would be required in order to produce powerful coronae), then (see also Taam & Lin, 1984; Burm, 
1985) 



-Pmag — QiOy PgasPtot, (2) 

where Qo is a constant (not necessarily smaller than unity, see M03). Then it is easy to show that, locally in the disc, the fraction of power 
released into the corona I Haardt & Maraschi 1991 Svensson & Zdziarski 1994 1 is given by 

f( R )~^(i+%*!&) 1/4 . (3 ) 



Pgas (R) 

The first important implication of the above result is that the coronal fraction / is reduced when radiation pressure dominates 1 . As radiation 
pressure is more and more dominant with increasing accretion rates, a dependency of the spectral properties of an accreting source (that are 
very sensitive to the fraction of power released in the hot, optically thin corona) on the accretion rate is introduced. This may be relevant, for 
example, for understanding spectral transition in galactic black hole binaries ( Merloni & Fabian 2002 1. 

At any given accretion rate the coronal emissivity profile F C (R), which determines the intensity of the hard X-ray radiation illuminating 
the disc at every radius, is given by 

F C (R) = f(R)Q(R) = ergs cm- s"\ (4) 

where Q, B and C are the relativistic correction factors of Novikov and Thorne (1973), approaching unity at large distances (see Appendix). 
One of the key feature of such a relativistic solution is the 'no-torque at the inner boundary' (in the following, NT) condition. Such a condition 
ensures that the disc radiative flux vanishes at the inner boundary, and this in turn implies that the ratio of radiation to gas pressure decreases 
as the inner boundary is approached from outside. In fact, the innermost part of a standard geometrically thin and optically thick disc is 
always gas pressure dominated if the NT condition is satisfied, and the coronal emissivity should therefore be strongly peaked towards the 
center. 

In order to show this, we have calculated the coronal emissivity profiles for a supermassive black hole (m = 10 s ) accreting at about three 
per cent of the Eddington rate (a rate such that at least a portion of the inner disc is dominated by radiation pressure) and with ao = 0.49, 
for different values of the black hole spin. We solve the relativistic disc structure equations, modified by the presence of the corona (see Eq. 
A2-A1 1), with the disc-corona coupling given by Eq. An example of the results is shown in Fig.^. The emissivity profile of the corona 
F c is much more strongly peaked than the corresponding disc emissivity profile Q. This effect does not depends on neglecting the effects of 
pressure gradients near the inner boundary. Indeed, as it is also clear from Fig.^, the peak of the coronal emissivity is not located near the 
innermost stable orbit, if "near" here means so close that substantial deviation from the values of the physical quantities computed with the 
Novikov-Thorne equations should be considered. In fact, the coronal emissivity profile peaks close to the radius where the gas pressure is 
equal to the radiation pressure. Calculations done for high accretion rate discs, including the effects of pressure gradients, show that the ratio 
of gas to total pressure increases substantially when going towards the center ( Abramowi cz et al. 198 8 I. from r/n n = 5/3 to r/rjn = 1 
(see e.g. Fig. 6 of Muchotrzeb & Paczynski, 1982, or Fig. 7 of Abramowicz et al. 1988), even for much higher accretion rates than those we 
are considering. Indeed, we can compute from our numerical solutions the ratio of the radial velocity at the location of the maximal coronal 
emissivity, v r (rmax) = u max , and compare it with the local sound speed, to have a measure of the distance between r m ax and the sonic 
point. Moreover, the error we make by neglecting pressure gradients can be estimated calculating the relative departure of the gas angular 
velocity from the Keplerian value at r = r max , from the momentum equation in the r direction: 

AQ (Q, 2 - ^k) 1/2 1 (l dP dv r , 

+ v r - f- (5) 



fix Ok ?"^k \p dr dr 

InFig.EJ) we show the location of the inner disc radius, of the radius of maximal coronal emissivity and of maximal disc emissivity 

1 Indeed, many of the results that are discussed here are valid in a more general framework than that proposed by M03. It is in fact sufficient that the magnetic 
pressure, Eq. |5J, be not proportional to the total pressure only in order for a dependency of / on the pressures ratio to be introduced. 
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Figure 2. The profiles of the coronal fraction f(R) for Ae = (NT, thick solid lines) and for the case of disc torque (DT) with Ae = 0.05, 0.1,0.5 (dashed 
lines, Ae increasing from top to bottom). The dimensionless black hole spin a = 0.99,0.9,0.7,0 is indicated in the different panels. The dimensionless 
accretion rate, defined as the ratio of the emerging (accretion) bolometric luminosity to the Eddington one has been kept fixed to rh = ^b.acc/^Edd = 
e(a)-M u c 2 /^'Edd = 0.03, where Mq is the constant external accretion rate and e(a) is the standard radiative efficiency computed in terms of binding energy 
at Rj n = Rms, which is a function of the dimensionless black hole spin a. 



without corona (upper panel) together with the values of u max /c s and AQ/£Ik (lower panel), as a function of the dimensionless black hole 
angular momentum a for the NT case. At the location of the maximal coronal emissivity, the radial flow is about three orders of magnitude 
subsonic and the error we make in neglecting pressure gradient and advective terms is a fraction of a percent. 

The coronal emissivity profiles shown in Figure^ represent an upper limit on the central concentration of hard illuminating X-rays that 
can be achieved without the need of any kind of energy extraction from matter inside the marginally stable orbit. Comparison of the calculated 
emissivity profiles with those needed to explain the broad relativistic lines observed in the spectra of the Seyfert 1 galaxy MCG-6-30-15 
(Wilms et al. 2001 . Fabian et al. 2002 1, reveals that the radial extent of the region with emissivity index (3 > 3 (the "core" region, see below) 
is probably too small in the NT scenario to account for the observations. In the next section we investigate the effects of applying a torque at 
the inner disc boundary on the coronal emissivity profiles. 



3 TORQUE AT THE INNER DISC BOUNDARY: THIN VS. THICK DISCS 

The original idea that magnetic field inside the radius of marginal stability (7? ms ) may alter the dynamics of accreting material and act to 
effectively torque the disc at its inner boundary was already mentioned in the early works on black hole accretion discs (e.g. in Page & Thorne, 
1974), although the standard hypothesis, in the following three decades of accretion disc research, has been that of no-torque at the inner 
boundary (see Abramowicz & Kato, 1989, and references therein). When the NT inner boundary condition is adopted, the maximal energy 
per unit mass available to be radiated by the accretion disc is uniquely determined by the binding energy at the innermost stable orbit and this 
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fixes the efficiency of accretion. On the other hand, if dynamically important magnetic fields in the plunging region exert significant stresses 
on the disc (disc torque condition, in the following DT), the available energy per unit mass increases and the efficiency may even reach values 
greater than unity, if magnetic stress can be exerted on the disc from material that has plunged into the ergosphere, thus effectively tapping 
the rotational energy of the black hole HCrolik 19991 lG ammie T999l|Agol & Krolik 2000) . As a consequence, the nature of the inner disc 
boundary condition, and in particular the extent of the torque exerted by the plunging gas on the accretion flow, depends on the complicated 
nature of the (turbulent) angular momentum transport in the disc and in particular on the level of the saturation magnetic field inside the 
disc (Reynolds & Armitage 200T]|Krolik & Hawley 2002) . The recent MHD simulations of Hawley and Krolik (2002), which are the ones 
with the best resolution to date, seem to confirm the idea that magnetized discs do not satisfy the no-torque at the inner boundary condition, 
showing that the stress generated by the MHD turbulence does not vanish near the marginally stable orbit. 

There is in fact a contrasting argument in such a debate. Indeed, all the simulations that have shown signs of magnetic stress inside the 
marginally stable orbit were non-radiative, and involved moderately thick discs (very thin, radiatively efficient discs are of course much harder 
to simulate). Thinner discs being denser, the Alfven speed is reduced in their inner part, and maintaining the causal connection between the 
magnetized gas in the plunging region and the disc should be more difficult. Indeed, Afshordi & Paczynski (2002), by studying a numerical 
model of a steady state, thin and isothermal accretion disc, argued that the NT inner boundary condition is appropriate, provided that the 
structure carrying the stress in the disc (likely the magnetic field) has a scaleheight much smaller that the radius, i.e. is itself geometrically 
thin (see also Li, 2003). 

It is worth stressing here that the very existence (and the variability properties) of broad and redshifted fluorescent Iron Ka lines in the 
spectra of accreting black holes, put constraints on the disc models: the discs that produce them have to be cold (not too highly ionized), 
optically thick to Compton scattering and illuminated by a strongly variable source of hard X-rays. Clearly, non-radiative thick accretion 
flows, as those simulated by Hawley and Krolik (2002), do not satisfy those conditions. 

Therefore, beside the no-torque (NT) and the disc-torque (DT) boundary conditions advocated, respectively, by Paczynski and collab- 
orators and by Krolik and collaborators, we envisage a third plausible alternative, applicable to the case in which a geometrically thin, cold 
disc is indeed present together with the hotter coronal phase, as inferred from the detection of broad iron Ka fluorescent lines. The proposed 
geometry for the inner flow is the following: a cold thin disc, with no torque at the inner boundary, giving rise to the observed broad lines 
via reflection of hard X-ray photons, is surrounded by a magnetically dominated (Merloni & Fabian 2001 1 corona. Magnetic field lines that 
penetrate the ergosphere at high latitude above the equatorial plane are connected to the coronal plasma, and the torque exerted by the rotation 
of the inertial frames transports energy into the corona (coronal thread 2 case, CT). This is in practice a generalization of the Blandford-Znajek 
(Blandford & Znajek 1977) mechanism, where the energy extracted from the rotating hole instead of powering a relativistic outflow is dis- 
sipated, via magnetic reconnection, in the corona. The coronal emissivity profile can then be more centrally peaked than in the NT case, 
because of the excess field dissipation in its innermost part. 

In the following sections we treat separately the two extreme case of a torque applied to the thin disc only (DT, section l3~"fl and of a 
direct magnetic connection between the black hole and the coronal plasma (CT, section l3~2l . First, we solve the equations for the coupled 
accretion disc-corona system and derive the coronal emissivity profiles in the two cases, then (section |4j we discuss their effects on the 
properties of the fluorescent iron lines profiles. 



3.1 Torque exerted on the thin disc 

We follow here the argument presented by Krolik (1999) and adopt the formalism of Agol & Krolik (2000), in particular, we assume nonzero 
stress at the marginally stable orbit, so that the flux at the disk surface in the fluid frame is not zero there. The relativistic equations for the 
disc structure discussed in Agol & Krolik (2000) are modified here by the inclusion of a term that takes into account the fraction of power 
that goes into the corona and supplemented with the coupling provided by Eq. The corona, in this case, is coupled only to the underlying 
disc and not to the free-falling gas, nor to the black hole horizon. 

Using this boundary condition, we have, for the locally generated radiation flux emerging from the disc: 

3 GMM (l-f) f xUsCUs^ 
8tt R? \ C(a;)x 1 /2 

Here Ae is the additional radiative efficiency relative to the one computed in terms of binding energy at r ms (the usual disc efficiency, 
e(a)), x = R/R g = Rc 2 /GM and C(x) is one of the auxiliary functions that incorporate relativistic corrections to disk equations 
(from Novikov and Thorne 1973, see also the Appendix). Rp T is a factor that represent the relativistic correction factors in the NT case 
(Novikov & Th orne 19731 ; it goes to zero at i? ms . The corresponding angular momentum conservation equation reads 

m q k (x) f x*£vl£i 

2tt I C{x)x 1 / 2 



F W = n>3 + R R • ( 6 ) 



T r<t> ( 2 ) dz = KZ I n,-S-1l* + T I ( 7 ) 



where in the above formula R^ T represent the torque relativistic correction factor (Novikov and Thorne, 1973), analogous to the one included 
in the flux formula. 



2 As our corona is not accreting, and no angular momentum conservation equation is discussed, it would be incorrect to talk of coronal torque. 
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Figure 3. The profiles of the coronal fraction f(R) for Ae = (NT, thick solid lines) and for the case of extra energy dissipation in the corona (CT) with 
Ae = 0.05, 0.1, 0.5 (dashed lines, Ae increasing from bottom to top). The dimensionless black hole spin a = 0.99, 0.7 and the index of coronal dissipation 
q e are indicated in each different panel. The dimensionless accretion rate, defined as the ratio of the emerging bolometric luminosity to the Eddington one has 
been kept fixed to rh = L rm f, i /L-^dd = e ( a ) Mqc 2 / L^dd = 0.03. 



With the help of Eqs.JSJ and 0, we can again extract the disc structure equations and find the solutions for the coronal fraction as a 
function of R, for different Ae. This is shown in Figure [2] for different black hole dimensionless spin parameters and for three different 
values of Ae, between 0.05 ( < e(0) ~ 0.06) and 0.5 ( > e(l) ~ 0.42). We have fixed the dimensionless accretion rate, defined as the ratio 
of the bolometric luminosity to the Eddington one m = I/b,acc/-i'Edd = t(a)Moc 2 / L^dd = 0.03. 

As expected, when a torque is applied to the geometrically thin, optically thick disc only, the excess dissipation at the inner boundary 
reduces the ratio P ga s/-P r ad and, consequently, /. The more energy is extracted from the plunging region and transferred to the disc, the less 
coronal emission is produced in the innermost region of the accretion flow. 

As it is well known I Agol & Krolik 2000] ILi 20 02 1, when extra torque is applied in the innermost region of an accretion disc, the 
emissivity profile of the disc itself 'becomes much steeper. Then, if the coronal fraction does not decrease too rapidly with decreasing distance 
from the hole, the coronal emissivity may be peaked, too. However, the intensity of the non-thermal coronal flux compared with the disc 
one is much reduced in the DT case, up to the point where coronal emission becomes negligible. In practice, given our prescription for the 
disc-corona coupling, a torque applied at the inner boundary of the thin disc only, has the effect of quenching the coronal activity in the 
innermost part of the accretion flow. In section |4| we will show a few quantitative measures of the amount of peakedness of the coronal 
emissivity profile, and discuss its observable effects on the iron Kq line profiles for different black hole spin parameters and Ae. Before that, 
we examine next the effect of extra energy deposition on the corona itself. 
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Figure 4. Plots of the ratio § of the core radius i?t, r to the radius of maximal coronal emissivity i? ma x (panel a), and of the ratio of the core coronal 
luminosity L CO rc (computed from iJ ma x to i?br)> to the outer one (from R^ T to Rout = 100-Rg = 200GM/c 2 ; panel b), both as functions of the black hole 
dimensionless angular momentum parameter a. The thick black solid lines represent values computed for Ae = (NT). Three disc torque cases are shown 
with (red) dashed lines for Ae = 0.05, 0.1 and 0.5, from bottom to top, respectively. CT cases correspond to (green) dot-dashed lines (q e = 4) and (blue) 
dotted lines (q e = 4.5), again going from Ae = 0.05, 0.1 to 0.5 from bottom to top. In all the calculations, we have fixed the dimensionless accretion rate to 
rh = 0.03. 



3.2 The black hole-corona connection 

The magnetically dominated plasma in the low-density corona (where Alfven speed may become relativistic) may indeed be the repository of 
most of the energy extracted by electromagnetic processes from a spinning black hole (Krolik 1999 Blandford 2001 Semenov et al. 2002 1. 
In this respect, we may consider the inner part of the magnetically dominated corona as the non-stationary analogous of the force-free 
magnetosphere envisaged by Blandford & Znajek (1977). The extra energy may also come from an increased coronal activity produced by 
the strongly magnetized plasma in the plunging region <Gammie* 1999 Krolik 1999 1, where efficient reconnection events may disconnect the 
free-falling plasma from the disc (thus reducing the effective torque at the inner disc boundary), while simultaneously increasing the energy 
available for coronal heating. The time-averaged effect of such a process of energy extraction and deposition (whose physical modeling, 
involving the electrodynamics of a reconnecting field in a curved background spacetime, is beyond the aim of this paper) on the radial profile 
of the hard X-ray radiation illuminating the underlying cold disc is the generation of a central cusp. 

In order to consider the effect of a black hole-corona connection, with the consequent increase of coronal emissivity in the simplest 
possible way, we apply a similar argument to that presented in the previous subsection. The effect of the excess energy dissipation in the 
coronal region, is parametrized by the same factor Ae representing the extra radiative efficiency of the accreting gas with respect to the 
standard value e(a). However, in this case, we leave unchanged the disc equation for the locally generated flux: 

3 GMMp(l — f*(R)) nt m 

as well as the angular momentum conservation equation: 

M n K (x) NT 

T r<t ,{z)dz = — R T , (9) 

which are therefore both identical to the no-torque case. 

Then, we modify the expression for the local coronal flux by simply adding a term with a power-law profile: 

P CT, n s SGMMp ( f*(R) QB u _ qc \ _ 2 _ x 

F c (R) = I —^3-3171 +K,eR I ergs cm s . (10) 

The additional term is normalized in such a way that the extra coronal luminosity L e = 2 2ivRQ e (R)dR, where we have defined 
Q e (r) = {3/%-K)GMM K e R~ qc , is equal to AeM c 2 , and thus we have 

/C e = ^(q c - 2)R^- 2 \Ae (11) 

Here we need to introduce a new definition of the fraction of power released into the corona. We define f(R) = F^ T (R)/Q(R) = 
f*(R) + f°(R), where now f*(R) is the (usual) fraction of dissipated accretion energy deposited into the corona (bound to be less than 
unity, see section|2j, that enters the disc structure equations, while f c (R) is the additional coronal energy in magnetic field lines connected 
to the plunging gas at moderate-to-high latitude. As opposed to the DT case, here the more the energy extracted from the black hole (i.e. 
the bigger Ae), the larger the coronal fraction f(R). Indeed, as expected, in the innermost disc region, the coronal fraction, i.e. the ratio 
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of the energy dissipated in the optically thin, hot, magnetically dominated corona to the total gravitational energy locally dissipated by the 
accretion flow, Q, may be larger than unity. This is clearly shown in Figure|5| where, again, we show the profiles of f(R) for both highly 
and moderately spinning black holes, a — 0.99 and a = 0.7, for two different values of the index q c and for three different values of Ae, 
between 0.05 and 0.5. Globally, we can define the integrated fractions of power released into the corona in the following way: 

^ 2f" n 2*Rr(R)Q(R)dR ^ 

> 2/~ 2nRQ{R)dR Lb,ac C ' 

Then we get, for the global value of the fraction of power released into the corona, < / >: 

<f>= U + Lo = *!+</•>, (13) 

^b.acc e(,1J 

from which we deduce that the extra term due to the energy transported from inside i? ms and released into the corona, < f c > is given by 
the ratio Ae/e(a) of the extra radiative efficiency to the standard (accretion) one. 



4 EFFECTS ON THE FE Ka LINE PROFILES 

In order to estimate the effect of the coronal emissivity profiles on the shape of the fluorescent Iron Ka line produced by the underlying cold, 
optically thick disc, we have approximated the computed F c (r) with a double power-law (see Fabian et al., 2002, for an application to the 
observation of MCG-6-30-15), assuming that the inner radius corresponds to the radius of maximal coronal emissivity (Ri n = -Rmax). For 
the inner core, extending from 7? m ax out to a break radius R^, we have approximated the peaked profile with an emissivity index j3 = 4, the 
same for all the computed profiles. Thus, for the normalized coronal flux we have 

Fc{R) = (^p) for <R<Rbr 

Fc{R) = ^K^y^y 2 for R bT < R < R out . (14) 

By fitting such a profile to the calculated F c , we find the core radius R^ T and the emissivity index, fc, of the outer, flatter part. 

Also, returning radiation I Cunningha m 1976) has to be taken into account 3 . This can be done approximately, by using the analytic fitting 
formula for the fraction of the returning radiation to the locally generated one of Agol & Krolik (2000). For R ou t ^ WRin such a fraction, 
R Te t can be approximated with a constant, that we denote i? rc t,oo- The returning flux then scales as r -3 : 

Fret (R ) ~ Rrct,co Q(Rout)(^) . (15) 



R 

We can then introduce the integrated luminosity from the core region L corc = 4-7T J^ hl [F C (R) + F rct (R)]RdR (and, correspondingly, that 

of the outer region L ou t = 4-7T J^ 011 ' [F C (R)) + F Ict (R)]RdR). By integrating these expressions, we get, for the relative luminosity of the 
two regions, (for fc 7^ 2) 

Lcorc _ fc - 2 ( (f3_l)+2A£V(£-l) 



(16) 



Lout 2 ^(l- x (2-«) + (p_2)^ 4 X 2 (x-l) / 

where £ = R hl /R ms .x, X = ftut/-Rbr and A = i? r ct,ooQ(-Rout). 

Both £ and L COIC /L ou t can be used as measures of the peakedness of the coronal emissivity profile. The values of these two quantities 
determines the profile of the fluorescent iron line produced by the reflection of the illuminating coronal radiation on the cold disc. In Figure^ 
we plot the ratio £ = Rhr/Rmax of the break radius to the inner radius as a function of the black hole dimensionless spin a, for different 
values of Ae, for the cases of DT, CT and NT. In Figure^ we plot the ratio Lcorc/Lout of the luminosity of the coronal core region to 
the luminosity of the outer coronal region. These figures show that the inner core in the DT case is less prominent than in the NT case for 
a < 0.99 and Ae < 0.1. Only for Ae>0.5 do the coronal emissivity profiles in the DT case become substantially peaked. On the other 
hand, for the coronal-thread case (CT), already for Ae = 0.05 is the ratio L COIC /L ou t larger than in the NT case. Its value increases rapidly 
for increasing Ae and q c . 

To compute the line profiles we have adopted the Laor (1991) model as implemented in XSPEC. Each line is the sum of two components, 
one produced by the core, the other by the outer region. The model for each component has the following parameters: inner and outer radius, 
emissivity index /3, line energy, inclination angle and overall normalization. In all cases, the line energy was fixed to 6.4 keV, the emissivity 
index of the core region to /3 — 4, the outer radius to 100i?s and the inclination angle to 30° . The variable parameters that determine the line 
shape are: the inner radius (coincident with i? ma x), the break radius, R^ r (which is simultaneously the outer radius of the core component 
and the inner radius of the outer component), the emissivity index of the outer component, fc, and the normalizations of the two components. 



3 In fact, also the vertical structure of the illuminating corona may become important, and a full general relativistic treatment of light bending effects 
1 Martocchia, Matt & Karas 2002 1 for realistic coronal geometries should be relevant to exactly determine the line shape in the most extreme cases. 
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Table 1. Physical parameters (inner radius, r max ; core radius, rt, r , both in units of Schwarzschild radii; outer emissivity index, @2', and relative normalizations 
of the core and outer component), obtained from a double power-law fitting of the coronal emissivity profiles for different values of the black hole dimensionless 
spin parameter a, of the extra radiative efficiency Ae, and of the inner boundary condition (T). In the fits, the core emissivity index has been fixed to 4 and the 
outer radius to 100-Rg- 

Model parameters Fit parameters 

T a Ae q e r max r br 2 N COIC N out 



NT 


0.99 







0.873 


0.897 


2.06 


0.068 


0.932 


NT 


0.9 







1.49 


1.78 


1.97 


0.169 


0.831 


NT 


0.7 







2.33 


2.88 


1.90 


0.194 


0.806 


NT 


0.1 







4.27 


5.40 


1.81 


0.219 


0.781 



DT 


0.99 


0.1 




0.734 


0.826 


2.11 


0.112 


0.888 


DT 


0.9 


0.1 




1.15 


1.25 


2.06 


0.086 


0.914 


DT 


0.7 


0.1 




1.67 


1.83 


2.01 


0.088 


0.912 


DT 


0.1 


0.1 




2.79 


3.02 


2.02 


0.077 


0.923 


DT 


0.99 


0.5 




0.734 


1.09 


2.30 


0.331 


0.669 


DT 


0.9 


0.5 




1.15 


1.51 


2.23 


0.241 


0.759 


DT 


0.7 


0.5 




1.67 


2.11 


2.25 


0.210 


0.790 


DT 


0.1 


0.5 




2.79 


3.22 


2.29 


0.138 


0.862 



CT 


0.99 





.1 


1 


0.734 


0.999 


2.57 


0.267 


0.733 


CT 


0.9 





.1 


4 


1.15 


1.62 


2.85 


0.293 


0.707 


CT 


0.7 





.1 


4 


1.67 


2.40 


3.02 


0.306 


0.694 


CT 


0.1 


0. 


.1 


4 


2.79 


4.30 


3.11 


0.361 


0.639 


CT 


0.99 


0. 


.5 


■I 


0.734 


2.34 


2.58 


0.689 


0.311 


CT 


0.9 


0. 


.5 


4 


1.15 


3.75 


2.74 


0.700 


0.300 


CT 


0.7 





.5 


4 


1.67 


5.06 


3.06 


0.679 


0.321 


CT 


0.1 


0. 


.5 


1 


2.79 


7.26 


3.44 


0.632 


0.368 


CT 


0.99 


0. 


.1 


1.5 


0.734 


1.24 


2.55 


0.411 


0.589 


CT 


0.9 





.1 


4.5 


1.15 


2.09 


2.88 


0.455 


0.545 


CT 


0.7 





.1 


4.5 


1.67 


3.20 


3.10 


0.484 


0.516 


CT 


0.1 


0. 


.1 


1.5 


2.79 


5.38 


3.31 


0.495 


0.505 


CT 


0.99 


0. 


.5 


4.5 


0.734 


3.28 


2.48 


0.780 


0.220 


CT 


0.9 





.5 


4.5 


1.15 


5.65 


2.58 


0.805 


0.195 


CT 


0.7 





.5 


4.5 


1.67 


8.79 


2.95 


0.823 


0.177 


CT 


0.1 





.5 


1.5 


2.79 


20.2 


3.49 


0.886 


0.114 



These have been fixed by imposing iVcorc + iVout = 1 and taking N core /N on t — Lcoro/iout as computed with the help of Eg. 1161 TableQ 
shows the value of the parameters of the double power-law fit to the coronal emissivity profiles that have been used as inputs into the Laor 
model used to produce the line profiles shown in figures|5||6|and0 

In Figure|5|we show the composite lines we obtain for the NT case for four values of the black hole spin (a = 0.1, 0.7, 0.9, 0.99). The 
red dashed profiles show the line components produced from the core region, while the black, dot-dashed ones those from the outer regions. 
Obviously, the higher the black hole spin, the more red-shifted is the core-component, which appear as a small excess at low energies. In 
Figure|5|we show the profiles derived in the CT case for two values of the parameter Ae. Comparison with the NT case demonstrate that, 
for Ae = 0.1, the redshifted core components are substantially weaker and more redshifted (because r max is smaller in the DT case, see 
Table0. For larger Ae, the broad red shoulder produced by the core of the coronal emission becomes more evident, especially for rapidly 
spinning black holes. Finally, in Figure we show the effect of the modified coronal emissivity profile in the case of black hole-corona 
connection (CT), for two different values of the power-law index of the radial profile of the extra energy dissipation, q e . As we expect from 
the increase of L COIC /L ou t with increasing Ae shown in Fig.|4jj, the strongly redshifted component of the line, associated with the coronal 
core is more and more prominent. 

The line profiles obtained in the CT with an extra efficiency of the order of 0. 1 (and in the DT case for much larger values of Ae), are 
indeed consistent with those observed by XMM Newton in the X-ray spectrum of the Seyfert 1 galaxy MCG-6-30- 1 5 l Fabian et al. 2002 1, 
indicating the possibility that dissipation of rotational energy of the black hole tapped by electromagnetic processes has indeed been observed. 
The calculation we have performed clearly demonstrate that the line profiles are more sensitive to energy of the gas in the plunging region 
(and/or in the ergosphere) being magnetically transported and dissipated into the corona rather than into the disc, in the sense that the same 
amount of extra luminosity (additional to the bolometric accretion one) produces much stronger effects in the CT case than in the DT one. 
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Figure 5. The line profiles calculated from Laor model in the no-torque at the inner boundary case (NT, Ae = 0), with two components: a core, between the 
inner radius and the core radius, with steep emissivity profile (J3 = 4, red dashed lines), and an outer region, from the core radius to the outer radius, with a 
more shallow emissivity profile. The solid lines represent the total line emission. Four different values of the black hole spin are shown. 




Figure 6. As in Figure|5] in the disc-torque at the inner boundary case (DT), for Ae = 0.1 (left) and 0.5 (right) and four values of the spin parameter a. 
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Figure 7. As in Figure[5] in the corona-torque at the inner boundary case (CT), for Ae = 0.1 (left) and 0.5 (right), q c = 4 (top) and q c = 4.5 (bottom), and 
four values of the spin parameter a. 



5 DISCUSSION 

The calculations we have presented in order to include the effects of different inner boundary conditions into observables are obviously 
oversimplified, and we discuss below the main uncertainties that may affect our conclusions. 

First of all, we only considered stationary solutions of the accretion equations. As clearly demonstrated by a number of simulations 
{Reynolds & Armitage 200T||Hawley & Kro lik 2002 1 the nature of the accretion flow at the inner boundary is intrinsically time dependent, 
and the angular momentum transport there, which is mediated by the filaments with the higher relative magnetic field I Reynol ds & Armitage 200 \\ , 
can vary strongly on the dynamical timescale. However the line profiles we aim at modeling are usually derived integrating the observed 
radiation over many dynamical timescales. In fact, a direct comparison of time-resolved spectra could shed some light on the complex time 
variability of the inner accretion flow. 

Also, we apply a (generalized) a stress model to describe angular momentum transport, which may be questionable if, as we assume, the 
MRI-driven MHD turbulence saturates due to buoyancy of magnetic flux tubes (Balbus & Papaloizou 1999 1. But our approach is essentially 
phenomenological, and the work presented here aims at understanding the most important pieces of information hidden in marginal effects 
on observable physical quantities (for example, in the iron line profiles) rather than solve a complex theoretical issue. In this respect, the 
advantage of the simplifying a prescription is to allow a relatively straightforward interpretation of the observations. 

Another reason of uncertainties is the identification of the inner edge of the reflecting cold disc (the "reflection edge", Krolik & Hawley 
2002) with the radius of the marginally stable orbit around the central black hole. As discussed in I Krolik & Hawley 2002 1, the true reflection 
edge, i.e. the place where the integrated surface density of the disc has declined enough that the disc becomes effectively transparent, may 
depend on the accretion rate (and on the fraction of power released into the corona), if magnetic stresses in the inner disc are such that the 
radial inflow accelerates smoothly inward rather that having a sharp transition at R ms - Numerical simulations of cold, thin, radiative discs at 
different accretion rates are clearly needed in order to better verify these hypothesis. 
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Finally, contrary to the cases of NT and DT inner boundary conditions, where, for a given small number of assumptions, a full self- 
consistent analytic model can be solved for the disc-corona system, the coronal-thread hypothesis is admittedly more speculative: we do 
not have a specific model for how the energy can be extracted from the black hole and deposited there (but see discussion at the beginning 
of section ing , and we do not take into account the problem of angular momentum transfer in the corona (our coronae are effectively non 
accreting). Once again, given all the theoretical uncertainties, we believe it is more appropriate to adopt the simplest assumption about how 
the energy extracted from the black hole energizes the corona (see section l3~2l and explore its observational consequences. 

Given all the above uncertainties, it is obvious that all three kinds of boundary condition we have adopted are highly idealized. None 
of them is likely to fully describe the physical condition in the innermost part of a turbulent disc surrounded by a magnetically dominated 
corona. Nonetheless, their predictions for the intensity of the coronal emission compared with the quasi-thermal disc one, for the coronal 
emissivity profiles and for the shape of the fluorescent Iron Ka lines are qualitatively very different in the three cases and relatively easy 
to distinguish. That is what makes the use of such idealized models worthwhile. In fact, the current observations provides us with enough 
information to constrain the models. The shape of the iron line is very sensitive to the coronal emissivity profiles. Both the spin of the black 
hole (through the amount of gravitational redshift the line photons suffer) and the extra efficiency Ae (through the relative intensity of the red 
wing of the line, which fixes the emissivity index and size of the core region) can in principle be measured. Simultaneous multi-wavelength 
observations can set limits on the relative fraction of coronal versus disc emission, therefore helping discriminate between the CT and DT 
cases. The full picture can then emerge after a careful comparison with existing models of all these pieces of informations. 

6 CONCLUSIONS 

Our work demonstrates that detailed spectral analysis of the most broadened and redshifted iron Ka lines from accreting black holes (and 
of their variability) is a fundamental tool with which to constrain both the geometry of the innermost part of the disc-corona system and the 
nature of the inner boundary condition of relativistic accretion flows. 

The emissivity profile to be considered in analysing the data, to which the line shape is very sensitive, must be that of the illuminating 
flux, and therefore of coronal radiation, rather than that of the disc. When the generation of a magnetically dominated corona is included 
self-consistently in the relativistic equations for standard accretion discs, the coronal emissivity profiles are usually strongly peaked in the 
inner part, as far as a portion of the disc is radiation pressure dominated. The extent of such an inner coronal "core" depends crucially on 
the inner boundary condition. Here we have explored with our calculation three different (idealized) prescriptions: a classical no-torque at 
the inner boundary condition; a magnetic connection to the plunging region of the cold, geometrically thin disc only (disc-torque) and a 
magnetic connection to the plunging region of the vertically extended, magnetically dominated corona only (coronal-thread). We show that, 
while in the NT case indeed the coronal emission is more centrally peaked than the disc one, this effect is much enhanced in the CT case, with 
the strongest deviation from the standard emissivity profile appearing in the innermost coronal region for larger black hole spin and larger 
amount of energy tapped from the plunging region (and/or the black hole itself). On the other hand, when the torque is applied directly to 
the thin disc (DT case), although the emissivity profiles may become again more centrally peaked (for quite large values of Ae), the coronal 
power is much reduced, and we expect the emission to be completely dominated by the quasi-thermal radiation coming from the optically 
thick disc. 

By approximating the exact coronal emissivity profiles with a broken power-law, we have shown how the detailed line profiles are 
sensitive to the location of the break radius and to the relative normalization of the inner and outer illuminating regions. These parameters 
depend strongly on the black hole spin, on the amount of extra energy deposited into the corona, and on the nature of the inner boundary 
condition itself. The recently observed extremely broad and skewed fluorescent lines in the X-ray spectra of the Seyfert 1 galaxy MCG-6- 
30-15 <Wilms et al. 200T1lFabianetar 2Q02 1, may therefore be read as an indication that the extraction of energy from spinning black holes 
is most likely mediated by the magnetic corona rather than by the geometrically thin, optically thick accretion disc. 
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APPENDIX A: RELATIVISTIC DISC STRUCTURE EQUATIONS WITH MODIFIED VISCOSITY LAW 

Following Novikov & Thorne (1973), we define 

A = l + a 2 x~ i + 2a 2 x~ 6 

B = l + ax~ 3 

C = 1 - 3x~ 2 + 2ax~ :i 

V = = 1 - 2x~ 2 + a 2 x~ 4 

c i i a 2 -4 , 2 -6 , o 4 -8 

£ = = 1 + 4a x — 4a x + 3a x 

Q = x(l 7 2ax-3)i/2 ~ Xo - L5fl 

3(xi — a) 2 (x — x\) 3{x2 — a) 2 (x — X2) 3(xz — a) 2 (x ~ x 3 ) 



xi(xi — xi)(xi - x 3 ) (xq-xi) x 2 (x 2 - xi) (x 2 - x 3 ) (xo-xi) x 3 (x 3 — xi)(x 3 - X2) {xo — x 3 ) 



(Al) 



where x = v/ Rc 2 /GM, xq = \J i? ms c 2 /GAf, and x\ = 2cos(l/3cos 1 (a) — t/3), x 2 = 2cos(l/3cos 1 (a) + 7r/3), X3 = 
— 2 cos(l/3 cos - (a)). 

(a) When magnetic stresses are proportional to the geometric mean of gas and total pressure I Merloni 2003 1, in the case of radiation 
pressure dominated discs, we obtain the following expressions for disc scaleheight, pressure, density (both calculated at the disc mid-plane), 
temperature and ratio of the radiation to gas pressure, respectively: 

h ~ 10.6m(l - f)A 2 B" i C 1/2 V- 1 £- 1 Q (A2) 
P ~ 3.15 x 10 18 ( a m)- 8/9 r- 8/3 rh« /9 ^- 16/9 5 8/9 C -4/9 £: 8/9 Q 8/9 dyne (A3) 
p ~ 1.29 x 10- 4 («n)- 8/9 r 3/9 m - 10/9 (l - ^-^-W^s/^ts/gg-io/g g cm -s (A4) 

T ~ 1.86 x 10 a (am)~ 2/9 r~ 2/3 ml /9 A~ 2/9 T>~ 1/9 £ 1/9 Q 2/9 K (A5) 
§^ ~ 7.97 x 10 5 (am) 2/9 r- 7/3 mJ 6/9 (l-/) 2 ^ 20/9 S- 4 O- 8/9 £:- 10/9 Q 16/9 , (A6) 

"gas 
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where the dimensionless parameters are defined as: r = R/Rs, rho = 0.06Moc 2 /L B dd and m = M/Mq. 



(b) In the gas pressure dominated part of the disc, instead we have: 

h ~ 2mxlO- 2 C 3/S \amr 1/W r 21/20 ml/^l-f) 1/w (A7) 

P ~ 1.88 x 10 18 e 3/a2 («m)- 9/10 r- 51/20 m^ /5 (l - f)' 1 ' 10 A' 1 B 9/5 V 2/5 £ V*Q-*t* dyne (A8) 

p ~ 20.3C 11/32 (am)- 7/1 ^- 33/2() (l-/)- 3/1 ^hf/ 5 ^- 1 S 3/5 O- 1/5 f 1/2 Q 2/5 gcm- 3 (A9) 

T ~ 5.6 x 10 8 r 1/4 («m)- 1/5 r- 9/10 m^ /5 (l - /) ^ B -^ v -i/5 Q 2/ 5 R (A1Q) 

§^ ~ 1.31 x 10 2 r 35/32 («"i) 1/10 r- 21/20 m^ /5 (l-/) 9/10 ^- 9/5 I)- 2/5 f- 1/2 Q 4/5 . (All) 

-fgas 



The value £ ~ 1 /3, the numerical that takes into account the difference in the vertical density profile that should enter in the definition of the 
disc scaleheight for a gas pressure dominated disc, can be found by imposing continuity of all physical quantities at the boundary between 
the two regions (see Merloni, 2001). 
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